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P.  EFACE 

A  cryogenic  storage  tank  on  a  space  mission  is  subjected  to  vari¬ 
ous  types  of  heating  such  as  solar  radiation,  planetary  albedo,  plane¬ 
tary  infrared,  and  spacecraft  radiation.  This  heating  leads  to  strati¬ 
fication  of  the  liquid  cryogen  and,  through  equilibrium  with  the  gas 
phase,  to  a  concomitant  increase  in  tank  pressure.  As  a  result  engi¬ 
neers  concerned  with  cryogen  management  on  space  missions  must  plan 
frequent  depressurization  cycles  (or  continuous  venting).  These  lead 
to  a  loss  of  pressurant  vapor,  and  in  a  very  low-gravity  environment 
can  result  in  losses  of  cryogen  by  entrainment.  To  reduce  the  number 
of  depressurization  cycles,  pumps  can  be  installed  in  the  storage  tank 
to  circulate  the  liquid  cryogen  and  minimize  vaporization.  In  some 
cases ,  mixing  the  cryogen  can  Increase  the  time-to-vent  so  that  the  mis¬ 
sion  can  be  completed  without  depressurization. 

This  Memorandum  considers  the  possibility  of  rotating  the  storage 
tank  to  mix  the  liquid  cryogen,  thus  decreasing  the  rate  of  tank  pres¬ 
sure  rise  and  increasing  the  tlme-to-vent.  This  method  of  mixing  has 
the  additional  advantage  of  having  no  moving  parts  within  the  tank. 

The  rotational  rates  required  to  initiate  mixing  are  low  and  fall  within 
typical  ACS  (attitude  control  system)  operating  regimes.  The  method 
also  provides  a  means  for  thermal  control  of  space-storables.  This 
Memorandum  is  part  of  continuing  RAND  studies  of  thermal  protection 
systems  for  space-storage  tanks  and  should  be  of  interest  to  military 
and  civilian  space  agencies  concerned  with  deep  space  missions. 

One  of  the  authors,  Ivan  Catton,  is  an  Assistant  Professor  In  the 
Department  of  Engineering,  University  of  California  Los  Angeles,  and  a 
Consultant  to  The  RAND  Corporation. 


SUMMARY 


A  method  is  presented  for  determining  the  onset  of  a  thermally  in¬ 
duced,  convective  mixing  motion  in  a  model  of  a  cylindrical  space-stor¬ 
age  tank  subjected  to  a  constant  heat  flux  at  its  outer  boundary.  The 
model  simulates  a  completely  filled,  liquid  cryogen  tank  with  shear- 
free  ends  rotating  in  a  low-gravity  environment.  The  governing  distur¬ 
bance  equations  reduce  to  a  self-adjoint  eigenvalue  problem  for  the 
critical  Rayleigh  number  (the  stability  criterion).  Using  the  Raylelgh- 
Rltz  technique  of  approximating  eigenvalues  with  an  equivalent  varia¬ 
tional  principle,  the  critical  Rayleigh  number  is  calculated  as  a  func¬ 
tion  of  the  Taylor  number  (the  ratio  of  Coriolis  forces  to  viscous  forces) 
and  the  aspect  ratio  (length-to-diameter  ratio)  of  the  cylindrical  tank. 
The  results  indicate  that  the  critical  Rayleigh  number  is  a  monotoni- 
cally  increasing  function  of  Taylor  and  a  monotonically  decreasing  func¬ 
tion  of  aspect  ratio. 

The  ‘simplifying  assumption  of  a  full  tank  with  shear-free  ends  is 
limiting,  but  it  does  not  change  the  qualitative  features  of  the  results. 
It  is  found  that  relatively  small  rotational  speeds  will  initiate  and 
maintain  a  thermally  destratifying  convective  motion  within  a  vessel 
filled  with  liquid  hydrogen.  Thus  rotational  destratification  can  be 
used  to  supplement  oi  to  substitute  for  the  more  conventional  mixing 
methods  proposed  for  space-storage  tanks.  The  degree  of  mixing  is  a 
function  of  the  difference  between  the  actual  rotation  rate  and  the 
initiation  rotation  rate.  The  present  method  of  analysis  cannot  be 
used  to  calculate  the  degree  of  mixing. 
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I .  INTRODUCTION 

Because  of  their  excellent  propulsive  characteristics,  cryogens 
are  used  as  primary  propellants  for  space  vehicles.  A  major  problem 
associated  with  unmixed  space-stored  cryogens  is  the  determination  of 
the  rate  of  tank-pressure  rise  due  to  external  heating  (e.g.  solar  rad¬ 
iation)  .  Typically  a  storage  tank  is  designed  with  some  upper  limit 
on  pressure  and  an  associated  tank -wall  thickness.  If  the  allowable 
pressure  is  Increased,  there  is  a  corresponding  increase  in  tank-wall 
thickness  and  tank  weight.  For  short  missions  the  tank  pressure  is 
maintained  by  continuous  vent.  However,  for  long-duration  missions 
the  loss  of  cryogen  through  continuous  vent  is  excessive. 

To  circumvent  this  problem  the  tank  pressure  is  allowed  to  build 
up,  and  its  rate  of  buildup  is  controlled  by  employing  appropriate  tank- 
wall  Insulation.  This  method  is  inefficient  since  in  a  zero-gravity 
environment ,  buoyant  forces  are  nonexistent  and  no  convective  mixing 
can  occur.  The  cryogen  becomes  thermally  stratified  and  the  maximum 
heat-storage  capability  of  the  cryogen  is  not  utilized.  At  present 
there  are  two  upper  limits  that  can  be  placed  on  the  rate  of  tank-pres¬ 
sure  rise.  A  maximum  upper  limit  is  obtained  if  it  is  assumed  that  all 
the  heat  entering  the  tank  goes  to  the  vapor.  If  it  is  assumed  that 

all  the  heat  entering  the  tank  goes  to  the  completely  mixed  liquid, 

* 

then  a  minimum  upper  limit  is  obtained.  Information  from  Holmes  was 
used  to  construct  the  two  tank-pressure  history  curves  represented  sche¬ 
matically  in  Fig,  1.  Since  the  thermal  stratification  problem  is  not 
completely  understood,  the  tank  designer  must  use  the  maximum  upper 
limit  curve  for  safety  reasons.  Thus  frequent  depressurization  cycles 
must  be  scheduled  to  Insure  that  a  space  mission  of  long  duration  will 
be  successful. 

If  the  cryogen  is  completely  mixed,  a  200-  to  300-day  mission  can 
be  accomplished  without  depressurizing  the  tank.  To  reduce  the  number 
of  depressurization  cycles,  pumps  can  be  installed  in  the  tank  to  cir¬ 
culate  the  liquid  cryogen  and  maintain  the  cryogen  at  a  uniform  temper¬ 
ature.  In  this  study  we  consider  the  possibility  of  rotating  the  tank 
to  initiate  and  maintain  a  thermally  destratifying  convective  motion 

* 

L.  A.  Holmes,  McDonnell**  Doug  las  Corporation,  Santa  Monica, 
California  (personal  communication) . 
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wlthln  the  vessel.  Rotational  destratlflcetlon  could  supplement  or 
substitute  for  the  more  conventional  methods  of  mixing. 

In  order  to  examine  the  proposed  method,  the  Important  features 
of  the  system  are  Incorporated  Into  a  simple  theoretical  model  of  an 
orbiting  propellant  tank.  A  fluid-filled  rotating  cylinder  (with  angu¬ 
lar  velocity  Q)  In  a  zero-gravity  environment  Is  subjected  to  a  con¬ 
stant  heat  flux  q  at  Its  outer  boundary.  A  convective  motion  will  be 
Initiated  within  the  cylinder  when  the  thermally  Induced  buoyant  forces 
(radially  directed)  are  sufficient  to  overcome  viscous  forces.  The 
marginally  unstable  state  Is  determined  by  the  critical  value  of  the 
Rayleigh  number.  Here  the  Rayleigh  number  Is  given  by 


^2,5 

R  .  (1) 

vkK 

where  b  Is  the  cylinder  radius,  and  a,  v,  k,  and  k  are  the  fluid  co¬ 
efficients  of  thermal  expansion,  kinematic  viscosity,  thermal  dlffuslv- 
Ity,  and  thermal  conductivity,  respectively.  The  critical  Rayleigh 
number  Is  a  function  of  a  parameter  which  appears  In  all  rotating  fluid 
systems ,  the  Taylor  number , 


(2) 


and  a  geometrical  shape  parameter,  the  cylinder  aspect  ratio  £/2b,  where 
£  Is  the  cylinder  length.  In  what  follows,  the  functional  relationship 
between  critical  Rayleigh  number,  Taylor  number,  and  cylinder  aspect 
ratio  Is  calculated  for  the  rotating  fluid  cylinder. 


o 
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II.  THE  GOVERNING  EQUATIONS 


Initially  a  quasi-incompressible  (Boussinesq)  fluid  fills  a  cir¬ 
cular  cylinder  which  is  rotating  about  its  z-axis  (see  Fig.  2)  in  a 
zero-gravity  environment.  A  constant  heat  flux  is  maintained  on  the 
cylinder's  radial  periphery  at  all  times.  The  initial  velocity,  tem¬ 
perature,  and  pressure  distributions  in  the  rotating  frame  of  refer¬ 
ence  are  given  by 


where  o  is  the  mean  fluid  density,  r  is  the  radial  coordinate,  and 
is  a  unit  radius  vector. 

It  is  assumed  that  instability  appears  via  a  marginal  stationary 
state;  hence  the  governing  perturbation  equations  in  dimensionless 
form  are 


div  u 

=  0 

curl  curl  u  +  grad  p  +  R 

rr^  +  \^  u  X  k  =  0 

(5) 

_2 

7  T  =  u 

•  re 

(6) 

where  u,  t,  and  p  are  the  velocity,  temperature,  and  pressure  distur- 

■”  2 
bances  measured  in  units  of  k/b,  qb/x,  and  ovk/b  ,  respectively.  Here 

k  is  the  unit  vector  along  the  z-axis.  If  it  is  assumed  that  all  per¬ 
turbations  are  axially  symmetric  (independent  of  the  azimuth  angle,  4), 
then  it  is  convenient  to  express  the  solenoidal  velocity  disturbance 
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as  the  superposition  of  a  poloidal  and  a  toroidal  vector  field.  Thus 


one  can  write 


u  =  curl  |^ijr(r,z)e^J  +  $(r,z)(^ 


where  if  and  $  are  defining  scalars,  and  ^  is  a  unit  vector  in  the  azi- 
muthal  direction.  After  eliminating  the  pressure  by  taking  the  curl 
of  Eq.  (5),  and  using  the  velocity  field  description  given  by  Eq.  (7), 
the  governing  perturbation  equations  assume  the  form 


curl  curl  curl  curl  (*e.)  +  R  curl  (rre  ) 

e  ~r 


-  T  curl  X  k)  =  0 


curl  curl  +  k  x  curl  =  0 


7  T  *  re  •  curl  ('^e  ) 

"HT 


Equations  (8) -(10)  may  be  written  in  an  equivalent  form  by  carry¬ 
ing  out  the  indicated  vector  operations.  One  then  obtains 


where 
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On  the  constant  heat  flux,  rigid,  radial  bounding  surface,  the  bound¬ 
ary  conditions  are 

^  °  ^ 

For  convenience,  the  cylinder's  end  bounding  surfaces  (z  =  ±  f/2b)  are 
assumed  to  be  Insulated,  shear-free  boundaries;  hence  the  boundary  con¬ 
ditions  on  these  surfaces  are 


i  =  ^=  M  =  |l=o 

^  ,2  dz  ^ 

dz 


at  z 


(15) 


Equations  (11) -(13)  and  the  boundary  conditions  given  by  Eqs.  (14) 
and  (15)  form  a  self-adjoint  eigenvalue  problem  for  the  Rayleigh  num¬ 
ber  with  the  Taylor  number  as  a  parameter.  A  variational  principle 
corresponding  to  the  eigenvalue  problem  can  be  constructed  as  follows: 
Multiply  Eqs.  (11),  (12),  and  (13)  by  i|f,  $,  and  t,  respectively,  and 
then  Integrate  over  the  cylindrical  region  (Indicated  by  {  )).  Inte¬ 
grating  by  parts  and  employing  the  boundary  conditions  yields 


=  0 


<».%)  -  (♦  i) 


Combining  these  expressions  yields  the  Rayleigh  number  as  the  extremal 
values  of 


= 


(16) 
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wlth 


L 


2 


9z 


(17) 


7 


2 


T  = 


(18) 


In  the  cylindrical  region  and  the  boundary  conditions  given  by  Eqs. 

(14)  and  (15).  The  Raylelgh-Rltz  approximation  method  Is  used  with 
the  above  variational  principle  to  obtain  upper-bound  estimates  to 
the  critical  Rayleigh  number.  A  trial  function  for  i|i  Is  constructed 
which  satisfies  all  the  boundary  conditions,  and  Includes  some  varia¬ 
tional  parameters.  Equations  (17)  and  (18)  are  then  solved  for  i  and 
T  subject  to  their  proper  boundary  conditions.  After  specifying  the 
Taylor  number,  the  Rayleigh  number  associated  with  the  trial  function 
can  be  found  from  Eq.  (16)  and  minimized  with  respect  to  the  parameters. 
The  accuracy  of  the  approximation  depends  on  the  nature  of  the  trial 
function  construction. 
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III .  THE  APPROXIMATE  SOLUTION  METHOD 


First  the  trial  function  sequence  i(f  is  constructed  from  a  com¬ 
plete  set  of  coordinate  functions,  each  of  which  satisfies  the  bound¬ 
ary  conditions  on  ijr.  It  is  convenient  to  choose  this  set  of  functions 
so  that  Eqs.  (17)  and  (18)  can  be  solved  with  little  effort  and  the 
integrals  appearing  in  Eq.  (16)  can  be  evaluated  analytically.  Thus 
\|l  is  chosen  as 


i|f  =  ^  Ajtj  (19) 

j=l 


where  A.  are  the  variational  parameters. 


j 


(r)  cos  az 


and 


J, (a.r)  I,  (a.r) 


F  ^ ,  J  -  .  ^  J 


Here  a  =  vb/i,  is  the  root  of 


Jj(a)l2(a)  +  J2(a)Ij(a)  =  0, 


(20) 


(21) 


J  is  the  Bessel  function  of  the  first  kind,  and  I  is  the  modified 
n  n 

Bessel  function  of  the  first  kind.  Using  the  above  expression  for  i|r 
and  Eq.  (17)  ,  one  can  obtain 


j  =  l 


sin  az 


(22) 


and 
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' '  A 


where 


2ot  I  (ar) 


a:z: . .  i  v?/’  i 


sin  az 


Substituting  ^  into  Eq.  (18)  yieldj 


_2  r 

V  T  =  ar  ^ 

j=l 


Equation  (25)  is  then  solved  for  t  using  the  method  of  variation  of 
parameters.  The  result  is 


■  I  V. 

j=l 


where 


Tj  -  a2|l__(.r)  r2f|(r)K^(ar)  a,  +  ^ 

-  K^(ar)  J  t^f^(r)I^(at)  drl  ,2,, 


and  Is  tha  modified  Bessel  function  of  the  second  kind.  Substitut¬ 
ing  the  expressions  for  *,  *,  and  t  into  E,.  (16).  and  minimising  with 
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respect  to  the  variational  parameters,  leads  to  the  following  infinite 
system  of  linear  homogeneous  equations  for  A.: 


(28) 


where  the  following  elements  of  the  coefficient  matrix  are  calculated 
in  the  Appendix: 


( 


0  0 

K)  ^  ^  -  2a  Qfj  - ^ ^  . 


\  4.4 

r  V  = 


j  =  k 


r»2  2  2 

4'  4 

■  “"j 


JoCOj) 

^  j  \ 


j  /  k 


/  2  ' 


,4  4a^aj  I^(a)  J^(a,) 


■‘--2-p-L  , 

^  r  X 


A  4 

aj  -  a^  Jj(aj) 


j  =  k 


,222 
=  4a  a, Of, 

J  k 


Io(a) 


O/x 


W 


TH  4\  TZ  4V^  4\  J,  (cyj 

-  a  j[a^  -  a  )  1  -  aj)(a  -  cy^j  1  j 


+  71;^ — 7xrz — 


j  k 


(29) 


The  element  Tj^  Tj^  cannot  be  expressed  in  any  simple  analytical  form; 
hence  it  is  computed  numerically.  This  calculated  element  and  the  above 
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two  elements  are  symmetric  with  respect  to  J  and  k.  The  characteris¬ 
tic  determinant  of  the  system  of  equations  is 


(30) 


For  a  given  Taylor  number  and  cylinder  aspect  ratio,  the  smallest  pos- 
sitive  root  of  Eq.  (30)  is  the  critical  Rayleigh  number  associated  with 
the  first  mode  of  instability.  If  the  trial  function  sequence,  Eq. 

(19),  is  truncated  after  a  finite  number  of  terms,  then  a  characteris¬ 
tic  determinant  with  a  finite  number  of  elements  is  obtained.  Increas¬ 
ing  the  number  of  coordinate  functions  in  the  approximating  sequence 
increases  the  accuracy  of  the  upper-bound  estimate  to  the  critical  Ray¬ 
leigh  number. 

Hie  following  criterion  is  used  to  decide  how  many  coordinate  func¬ 
tions  should  be  included  in  the  trial  function  sequence.  For  a  given 
Taylor  number  and  cylinder  aspect  ratio,  if  the  critical  Rayleigh  num¬ 
ber  associated  with  the  (n  +  l)-term  sequence  is  less  than  one  percent 
smaller  than  the  critical  Rayleigh  number  associated  with  the  n-term 

2 

sequence,  then  the  characteristic  determinant  is  truncated  with  (n  +  1) 
elements.  It  is  found  for  cylinder  aspect  ratios  of  interest, 

0.25  5  Jt/2b  5  4.0,  and  for  Taylor  numbers  less  than  10^^,  a  ten-term 
sequence  always  satisfies  the  above  criterion.  Thus  the  first  ten  roots 
of  Eq.  (21)  are  needed  for  the  calculations.  They  are  presented  in 
Table  1. 

The  computed  critical  Rayleigh  number  approximations  for  various 
Taylor  numbers  and  cylinder  aspect  ratios  are  presented  in  Table  2. 

These  results  are  graphically  illustrated  in  Fig.  3.  It  can  be  seen 
that  for  large  values  of  the  Taylor  number,  the  critical  Rayleigh  num¬ 
ber  is  a  linear  function  of  the  Taylor  number.  One  observes  that  R^/T 
is  a  monotonically  increasing  function  of  the  Taylor  number  and  ap¬ 
proaches  an  asymptotic  limit.  From  Eqs.  (1)  and  (2)  one  has 


R 

c 

T 


(31) 


Taylorr 

Number 


Table  2 

CRITICAL  RAYLEIGH  NUMBER  APPROXIMATIONS 


Cylinder  Aspect  Ratio 


1.2030 

X 

10 

1.2052 

X 

i 

10 

1.2265 

X 

1 

10 

1.4360 

X 

1 

10 

3.2359 

X 

lo' 

1.4833 

X 

10' 

8.5886 

X 

lo' 

5.5396 

X 

10* 

4.3031 

X 

10 

4.0501 

X 

10* 

4.0324 

X 

lo' 

0.50 

1.00 

2.00 

3.9241  X  10^ 

2.6329  X  10® 

2.3650  X  10-^ 

3.9391  X  10^ 

2.6388  X  10® 

2.3667  X  10® 

4.0886  X  10^ 

2.6972  X  10® 

2.3835  X  10® 

5.5512  X  10^ 

3.2752  X  10® 

2.5509  X  10® 

1.7766  X  lo'^ 

8.5341  X  10® 

4.1731  X  10® 

8.8904  X  lo'^ 

4.2656  X  lo'*^ 

o 

X 

5.0319  X  10^ 

2.3878  X  10^ 

9.2694  X  10 

3.1047  X  10^ 

1.4478  X  10® 

5.5015  X  10® 

2.1754  X  10^ 

9.5190  X  10® 

3,4873  X  10® 

1.8882  X  10® 

7.3645  X  10^ 

2.4408  X  10^ 

1.7733  X  10^ 

6.8208  X  10® 

2.0714  X  10® 

2.3018 

X 

10 

2.3022 

X 

10 

2.3066 

X 

10 

2.3501 

X 

10 

2.7820 

X 

10 

6.7770 

X 

10 

3.3984 

X 

10 

1.9281 

X 

10 

1.1879 

X 

10 

7  .7683 

X 

10 

5.8368 

X 

10 

relation  between  critical  Rayleigh  and  Taylor  numbers  for  various  aspect  ratios 
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Slnce  cyv/4kK  is  a  function  of  fluid  properties  only,  there  Is  a  limit 
on  the  heat  flux-radius  product  below  which  no  convective  flow  can  oc¬ 
cur  at  any  rate  of  rotation.  Figure  4  shows  this  lower  limit  as  a 
function  of  storage  tank  dimensions  and  cryogen  properties. 


Fig. 4— The  lower  limit  on  the  heat  flux-radius  product 
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IV.  DISCUSSION  OF  RESULTS 


Several  approximations  have  been  incorporated  into  the  model  used 
to  examine  the  feasibility  of  rotating  a  cryogenic  storage  tank  to  ini¬ 
tiate  a  convective  mixing  motion  within  the  liquid  cryogen.  The  cylin¬ 
drical  vessel  is  assumed  to  be  completely  filled.  This  is  almost  never 
true  (except  for  supercritical  storage)  since  an  ullage  must  be  avail¬ 
able  to  permit  thermal  expansion  of  the  liquid.  When  a  partially  filled 
vessel  is  rotating,  the  fluid  will  arrange  itself  so  that  the  ullage 
will  be  along  the  axis  of  rotation.  The  physical  problem  is  really  that 
of  fluid  contained  between  rotating  concentric  cylinders.  The  outer 
cylinder  is  rigid  and  the  inner  cylinder  is  a  shear-free  boundary.  Tnis 
problem  is  not  more  difficult  in  principle;  however,  the  mechanics  of 
solving  it  are  more  involved.  In  the  similar  problem  of  the  stability 
of  couette  flow,^^^  the  stability  criterion  increases  as  the  annular 
gap  width  decreases.  The  same  qualitative  result  is  expected  in  the 
rotational  destratification  problem. 

A  second  approximation  incorporated  into  the  model  is  to  assume 

the  cylindrical  tank  end  walls  to  be  shear-free.  This  approximation 

(2) 

is  made  by  Catton  and  Edwards  for  a  somewhat  different  physical  prob¬ 
lem,  but  one  with  essentially  similar  mathematical  characteristics.  The 
results  indicate  that,  if  the  end  walls  are  rigid,  the  additional  vis¬ 
cous  dissipation  that  must  be  overcome  is  inconseq^^antial  when  the 
tank  aspect  ratio  is  greater  than  2.0.  At  a  tank  aspect  ratio  of 
the  theoretically  predicted  stability  criterion  is  10  percent  lowe 
the  experimentally  observed  value.  Even  at  low  aspect  ratios  the  quali¬ 
tative  features  of  the  solution  are  valid  and  the  feasibility  of  ini¬ 
tiating  a  convective  flow  remains. 

Since  the  possibility  of  Inducing  a  convective  flow  in  a  heated, 
spinning  cylinder  has  been  demonstrated  theoretically,  it  is  useful  to 
consider  a  specific  cryogen  that  is  used  in  space  vehicles,  i.e.  liquid 
hydrogen.  The  thermodynamic  properties  of  liquid  hydrogen  at  37°R  given 

(3) 

by  Johnson  are  presented  in  Table  3.  Using  these  properties  and  Fig. 
3,  a  relation  among  heat  flux,  rate  of  rotation,  and  aspect  ratio  for 
a  liquid  hydrogen  tank  can  be  developed  as  shown  in  Fig.  5. 


Property 


Value 


As  an  example  of  the  application  of  these  results,  consider  a 

super-insulated  liquid  hydrogen  tank  with  a  heat  flux  of 
-1  -2 

0,1  BTU  hr  ft  ,  a  15-ft  radius,  and  a  60-ft  length.  It  is  found 
from  Fig.  5  that  a  rotational  rate  of  lO”^  sec”^  (one  revolution  every 
175  hours)  will  initiate  a  convective  flow  within  the  tank.  This  rate 
of  rotation  is  extremely  low  and  should  be  easy  to  achieve.  For  typi¬ 
cal  heat  fluxes,  tank  dimensions,  and  cryogens ,  the  spin  rates  required 
to  initiate  rotational  destratification  rail  within  typical  ACS  (atti¬ 
tude  control  system)  operating  regimes. 

This  Memorandum  has  considered  only  the  feasibility  of  inducing 
a  mixing  motion.  The  tank  would  have  to  be  rotated  faster  than  its 
critical  speed  to  achieve  a  specified  mixing  rate.  Experience  with 
similar  problems  indicates  that  increasing  the  rate  of  rotation  one 
order  of  magnitude  above  the  critical  speed  will  increase  the  effective 
thermal  conductivity  of  the  fluid  by  a  factor  of  five.  Further  increase 
in  spin  rates  will  increase  the  effective  thermal  conductivity  by  ap¬ 
proximately  the  square  root  of  the  spin  rate.  Since  the  ability  of  the 
fluid  to  transfer  momentum  (i.e.,  to  mix)  is  equivalent  to  the  ability 
of  the  fluid  to  transfer  heat,  mixing  rates  should  depend  on  spin  rates 
in  a  similar  manner.  It  is  not  known  at  this  point  how  to  put  a  fig¬ 
ure  of  merit  on  mixing  rate.  Even  if  it  is  required  that  the  spin  rate 


Fig. 5—Rotatlon -heat  flux  relation  in  liquid  hydrogen  for  various  aspect  ratios 
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be  several  orders  of  magnitude  above  the  critical  spin  rate  to  achieve 
adequate  mixing,  the  spin  rates  will  only  be  of  the  order  of  one  revo¬ 
lution  per  day. 

Examination  of  Fig.  5  reveals  an  interesting  paradox.  It  can  be 

calculated  for  the  above  example  that  if  the  heat  flux  is  reduced  to 

-1  -2 

less  than  0.035  BTU  hr  ft  ,  convective  mixing  will  not  occur  at  any 
spin  rate  because  the  centrifugal  body  forces  cannot  overcome  the  vis¬ 
cous  and  Coriolis  forces.  The  paradox  suggests  that  it  might  be  pos¬ 
sible  to  optimize  a  particular  mission  by  removing  tank -wall  insulation 
in  order  to  ensure  adequate  mixing  and  thermal  destratification.  The 
average  temperature  of  the  liquld-^'apor  mixture  would  rise,  but  the 
tank  vapor  pressure  would  fall.  A  study  of  the  trade-offs  among  insu¬ 
lation  thickness  (weight),  heat  flux,  boil-off  liquid  cryogen  weight, 
ACS  weight,  etc.  for  various  missions  would  be  useful  to  determine  the 
optimum  set  of  conditions. 

There  is  a  need  to  study  ways  of  starting,  maintaining,  and  stop¬ 
ping  rotation  of  the  large  tanks  used  in  long-duration  space  missions. 
The  example  tank  is  filled  with  almost  200,000  pounds  of  liquid  hydro¬ 
gen.  Add  to  this  an  equivalent  tank  weight,  and  it  is  clear  that  even 
the  maintenance  of  small  spin  rates  would  require  large  amounts  of  en¬ 
ergy.  The  start-up  problem  requires  a  study  of  the  time  it  takes  the 
cryogen  to  spln-up  to  the  tank's  design  speed.  Since  the  kinematic 
viscosity  of  cryoy,ens  is  low,  the  spin-up  time  may  be  long  and  the  tank 
designer  may  need  to  put  baffles  inside  the  tank  to  achieve  rapid  spln- 
up.  Once  the  tank  designer  decides  on  an  operating  spin  rate,  he  must 
examine  how  this  spin-rate  influences  the  ullage  and  subsequent  slosh¬ 
ing  and  surface  tension  effects.  These  are  nonlinear  problems  that 
are  beyond  the  scope  of  this  feasibility  study. 
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